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ABSTRACT 

Many supernova remnants (SNRs) are characterized by a knotty ejecta structure. The Vela 
SNR is an excellent example of remnant in which detached clumps of ejecta are visible as 
X-ray emitting bullets that have been observed and studied in great detail. We aim at mod- 
elling the evolution of ejecta shrapnel in the Vela SNR, investigating the role of their initial 
parameters (position and density) and addressing the effects of thermal conduction and radia- 
tive losses. We performed a set of 2-D hydrodynamic simulations describing the evolution of 
a density inhomogeneity in the ejecta profile. We explored diff'erent initial setups. We found 
that the final position of the shrapnel is very sensitive to its initial position within the ejecta, 
while the dependence on the initial density contrast is weaker. Our model also shows that 
moderately overdense knots can reproduce the detached features observed in the Vela SNR. 
Efl&cient thermal conduction produces detectable eff'ects by determining an efficient mixing 
of the ejecta knot with the surrounding medium and shaping a characteristic elongated mor- 
phology in the clump. 

Key words: Hydrodynamics - Shock waves - Methods: numerical - ISM: supernova rem- 
nants - ISM: kinematics and dynamics - ISM: individual object: Vela SNR 



1 INTRODUCTION 



The ejecta in supernova remnants (SNRs) drive the exchange of 
mass and the chemical evolution of the galactic medium. The 
structure of SNR ejecta has been proved to be knotty, and several 
, clumps have been observed at different wa yelegths in remn ants of 
■ core-collapse supemovae, as G292.0+L8 J Park et al .112004 1). Pup- 
pis A Ikatsuda et al. 2008), and Gas A, wher e knots have been 
observed also beyond the main shock fron t (iFesen et al.l l2006l 
iHammell & Fesenil2008l . [DeLanev et al.ll201Ql) . 

The Vela SNR, being the nearest SNR, represents a privileged 
target for this kind of studies, since it is possible to observe fine 
structures down to small physical scales. Despite the bulk of the 
X-ray emission of the Vela SNR is associated with the post-shock 
interstellar medium. X-ray emitting ejecta have also been observed. 
In particular, six protruding features, with characteristic boomerang 
morphology, (labelled Shrapnel A-F) have been identified in the 
Vela SNR by Aschenbach et al. (1995), who argued an ejecta ori- 
gin for these structures which appear to be detached from the rem- 
nant. The association with ejecta fragments has been supported by 
more recent observations performed with Chandra, XMM-Newton, 
and Suzaku. The analysis of the XMM-Newton observation of 
Shrapnel D, has revealed th at O, Ne, and Mg abundan ces are sig- 
nificantly larger than solar (iKatsuda & Tsunemill2005h . A similar 
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abundance pattern has been observed with S uzaku in Shrapnel B 
jYamaguchi & Katsudall 20091) . but in this case the overabundances 
of the lighter elements are less prominent, suggesting more effec- 
tive mixing with the interstellar medium (ISM). A Chandra obser- 
vation of Shrapnel A, whose projected distance from the center of 
the remnant is larger by - 20% than Sh rapnel D, reveals instead 
oversolar Si:0 ratios dMivata et al.ll200ll). Significant Si overabun - 
dance (Si~ 3) has been confirmed by iKatsuda & Tsunemil (|2006|) . 
who analyzed an XMM-Newton observation, finding solar or sub- 
solar values for the O, Ne, Mg, and Fe abundances. These results 
show differences in the chemical composition between Shrapnel 
A and Shrapnel B and D. In the northern rim of the Vela shell, 
iMiceli etalTfeOOS !) discovered new X-ray emitting clumps of ejecta 
whose projected position is behind the main shock front. The rel- 
ative abundances (0:Ne:Mg:Fe) of these new shrapnel are in good 
agreement with those observed in Sh rapnel D. Simil a r abun dance 
pattern have been observed also by iLaMassa et al.l (l2008h . who 
found ejecta-rich plasma in the direction of the Vela X pulsar wind- 
nebula. 

The present day morphology of SNRs and the structure of 
ejecta are believed to reflect the physical characteristics of the SN 
explosion (e. g., intrinsic asymmetries of the explosion, interac- 
tion of the early blast with the inhomogeneities of the circumstellar 
medium, physical processes in the aftermath of the explosion, etc.) 
and their detailed study promises to contribute to our understanding 
of the SN explosion physics. In the light of these considerations, it 
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is then interesting to model the evolution of the ejecta knots to un- 
derstand how the current position and chemical properties of the 
shrapnel in the Vela SNR depend on the physical conditions at the 
supernova explosion and on the dynamics of the explosion itself. 

The evolution of dense, supersonic clumps o f SN ejecta run- 
ning i n a uniform medium has been studied by I Anderson et alJ 
(1 1994*) and I Jones et alJ (Il994l) . who identified three main stages of 
evolution: a bow-shock phase, an instability phase and a dispersal 
phase. However, these models do not describe in detail the inter- 
action of the clump with the remnant (post- shock medium, main 
shock, reverse shock) and do not include impo rtant physical eff'ects 
(as thermal conduction and radiative cooling). [Cid-Femandes et all 
([1996) included radiative losses in their 2-D models, but they fo- 
cus sed on the interaction of a knot with a very small supernova 
remnant (~ 10^^ cm, i. e. more than 100 times smaller than the 
Vela SNR) evolving in an extremely dense medium (10^ cm~^). 
A hydrodynamic model (without thermal conduction and radiative 
coolin g) specifical ly tuned for the Vela SNR has been developed 
bv lWang & Chevaliej (120021) (hereafter WC02) who followed the 
evolution of a shrapnel by using 2-D simulations in spherical coor- 
dinates (because of the geometry of their simulations, the shrapnel 
are modeled as toroidal structures with very large masses). WC02 
did not model the early evolution of the ejecta knot, but started their 
simulations at the time trev, corresponding to the first interaction 
of the shrapnel with the reverse shock front. They explored diff'er- 
ent values of tyev and of the density contrast between the shrapnel 
and the surrounding ejecta, and found that, in order to produce 
an observable protrusion on the shock front (like that observed in 
Shrapnel A-F), a very high density contrast (x ~ 1000) is nec- 
essary. With lower density contrasts ix < 100), the shrapnel are 
rapidly decelerated and fragmented by hydrodynamic instabilities 
and the observed features cannot be reproduced (for theeffects_of 
hydro dyn amic instabilities on shocked clouds see also iKlein et aP 
1 19941 and lOrlando et al.|[2005h . Large density inhomogeneities in 
the clumps are difficult to explain in a core-collapse SN explosion. 
WC02 argued that a model that includes the eff'ects of radiative 
cooling may show that lower values of x are needed to match the 
observed protrusions. Also, WC02 do not include in their model 
the eff'ects of thermal conduction that, as shown by lOrlando et al.l 
12005), can efficiently suppress the hydrodynamic instabilities, thus 
allowing the shrapnel to overcome the main shock-front without be- 
ing disrupted. Recently, the evolution of knotty e jecta in a Type la 
SNR has been modelled bv lOrlando eTaP (l2012h (hereafter 012), 
who found that small clumps with initial x ^ ^ can reach the SNR 
shock front after ~ 1000 yr. Nevertheless, these ejecta knots are 
then rapidly eroded and do not produce significant protrusions in 
the SNR shock front, thus being unable to reproduce the features 
observed in the Vela SNR. 

Here we present a set of 2-D hydrodynamic simulations of 
the evolution of an (initially spherical) ejecta shrapnel in the Vela 
SNR. We include in our model both thermal conduction and ra- 
diative cooling and explore diff'erent values of x and of the initial 
position of the shrapnel in the ejecta profile. We aim at addressing 
the role of thermal conduction and radiative cooling and at under- 
standing how the initial properties of the shrapnel influence its evo- 
lution. We also aim at evaluating whether values of x lower than 
100 = 10 - 50) can reproduce the observed features. 

The paper is organized as follows: the hydrodynamic model 
is described in Sect. |2l the results of the simulations are shown in 
Sect. [3 and our conclusions are discussed in Sect.|4] 



2 HYDRODYNAMIC MODELING 

2.1 Initial conditions and model equations 

We model the evolution of a shrapnel in a SNR by performing a 
set of 2-D simulations in a cylindrical coordinate system (r, z), 
assuming axial symmetry. The system setup consists of a spher- 
ically symmetric distribution of ejecta with initial kinetic energy 
K - 10^^ erg (the initial thermal energy is only the 0.2% of K) and 
mass Mej = 12 Mq (representing the initial blast wave), where we 
place a dense, spherical, knot (the shrapnel) in pressure equilibrium 
with the surrounding ejecta and with central coordinates (0, Rs). 

The radial density profile of the ejecta consists of two power- 
law segments (p oc r~'" on the inside and p oc on the outside), 
in agreem ent with th e densi ty structure in a core-collapse SN de- 
scribed by IChevalied (l2005h . For our simulations, we use m = 1 
and b = 1 1.2, and the position of the transition between the flat and 
steep regimes is derived by following Chevalier (2005). The initial 
velocity of the ejecta increases linearly (up to 6 x 10^ cm/s) with 
their distance from the center. The maximum velocity is reached 
at the initial radius of the ejecta, i. e., = 4.5 x 10^^ cm. These 
values correspond to ~ 240 yr after the explosion, appropriate for 
the relatively late stages of the SNR evolution that we address in 
this study. In fact, the starting time of our simulations corresponds 
to only ~ 2% of the Vela SNR age and the shrapnel reaches the 
reverse shock ~ 2500 - 3000 yr after the explosion in all our simu- 
lations (i. e., the system has enough time to evolve, before the inter- 
action of the shrapnel with the SNR reverse shock occurs). We then 
conclude that our simulations can provide a realistic description of 
the actual conditions in the Vela SNR. 

The initial mass of the shrapnel is 0.05 Mej (1.19 X 10^^ g), 
its density is x times larger than that of the surrounding ejecta at 
distance Rs from the centeiQ, and its velocity is the same as that of 
the ejecta at distance Rg from the center. We aim at showing that 
the detached shrapnel observed in the Vela SNR can be the result 
of moderately overdense clumps of ejecta originating in relatively 
internal layers. We explored different values of ;^ and of Rs, namely 
X = 10, 20, 50, and Rs = 1/6 R^., 1/3 R^.. Figure [T] shows the 
initial density and temperature conditions for the case with;^ = 20 
and Rs - 11?) R^^j. The simulation setups discussed in this paper are 
summarized in Table [T] 

Vela SNR is the result of a core-collapse SN and we expect 
the ambient medium to be "perturbed" by the wind residuals of 
the massive progenitor star. However, we assume for simplicity a 
uniform ambient medium as in WC02, since here we are not in- 
terested in modeling the details of the remnant evolution. The fina l 
(i. e. after 11000 yr, the age of the Vela SNR. iTavlor et al.lll993h 
radius of the remnant strongly depends on the choice of the ambi- 
ent density value, njsM- We set njsM - 0.5 cm~^, because with this 
value (and with the chosen values of K, Mej, njsM, b, and m), the 
radius of the shell after 11000 yr is Rsheii ~ 5.4 x 10^^ cm, in good 
agreement with the observed radius of the Vela SNR, that ranges 
between 5x10^^ cm and 5.5 x 10^^ c m (by assuming a distance of 
250 pc, in agreement with Bocchino et al. 1999, Cha et al. 199^). 

Our model solves the time-dependent compressible fluid equa- 
tions of mass, momentum, and energy conservation. In three cases 
we ran the same simulation with/without thermal conduction and 
radiative cooling inside our system, as shown in Table [T] As for 
thermal conduction, we considered both the Spitzer and the satu- 
rated regimes, while radiative losses (that can play an important 
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Figure 1. Density (left panel) and temperature (right panel) 2-D cross- 
sections through the (r, z) plane showing an example of the initial condi- 
tions of our simulations (namely, Rl/3 - CHI20, see Table [T). The system 
consists of an expanding spherically symmetric distribution of ejecta where 
we place a dense, isobaric, and spherical knot. In this case the knot is 20 
times denser than the surrounding ejecta. 

role in the Vela SNR, as shown by iMiceli et"aDl2006l) were com- 
puted for an optically thin therma l plasma. The model equations are 
described in iMiceli et aP ([2006') (equa tions 1-5 therein) and were 
solved by using the FLASH code (Frvx ell et aPE oOQ). 

The computational domain extends over 8x10^^ cm in the r 
and z directions. We use axisymmetric boundary conditions at r = 0, 
reflection boundary conditions at z = 0, and zero-gradient (outflow) 
boundary conditions (for v, p, and p) elsewhere. We trace the mo- 
tion of the ejecta material and of the shrapnel with passive tracer^ 
Considering the large range in spatial scales of our simulations, 
we exploited the adaptive mesh capabilities of the FLASH code 
by adopting up to 10 nested levels of resolution (the resolution in- 
creases by a factor of 2 at each level). The refinement /derefinement 
criterion (Lohner 1987) follows the gradients of density, tempera- 
ture, and tracers. The finest spatial resolution is 1.95 x 10^^ cm at 
the beginning of the simulation, therefore there are 230 computa- 
tional cells per initial radius of the ejecta, and ~ 10 cells per initial 
radius of the shrapnel (that varies in the range 1 .8 - 2.9 x 10^^ cm). 
Because of the expansion of the system, the resolution is reduced 
by a factor of 2 after 2500 yr. We verified that by changing the res- 
olution of our simulations by a factor of 2, the results do not change 
significantly (see Appendix lAl for further details). 



3 RESULTS 

3.1 Evolution of the system 

We first focus on simulation RI/3-CHI20. Figure|2]shows the 2-D 
cross-sections through the (r, z) plane of temperature and density at 
diff'erent evolutionary stages of the 1/3 - CHI20 simulation. The 
left panel shows the system 5000 yr after the beginning of the sim- 
ulatior0, when the shrapnel interacts with the inter- shock region. 

^ Both tracers have zero mass and do not modify the dynamics of the sys- 
tem. 

^ In the following, all the ages will be reckoned from the beginning of our 
simulations. We remind the reader that our initial condition corresponds to 
~ 240 yr after the SN explosion. 



Table 1. Physical parameters of the model setups (see Sect.|2)- TC-RL in- 
dicates that the simulations were perfomed by including thermal conduc- 
tion and radiative losses. In all the setups the initial radius of the ejecta is 
R^. = 4.5x10^^ cm. 



Model setup 




X 


Rs (cm) 


TC-RL 


Rl/3-CHnO 




10 


1.5 X 10^^ 


No 


Rl/3-CHm- 


TR 


10 


1.5 x lO^^ 


Yes 


R1/3-CHI20 




20 


1.5 x 10^^ 


No 


R1/3-CHI20- 


TR 


20 


1.5 X 1018 


Yes 


RI/3-CHI50 




50 


1.5 X 10^8 


No 


R1/6-CHI50 




50 


0.75 X 10^8 


No 


RI/6-CHI50- 


TR 


50 


0.75 X 10^8 


Yes 




-26 -24.75 -23.5 -22.25 -21 

Log p [gm cm"^] 



Figure 2. Left panel: density (left) and temperature (right) 2-D cross- 
sections through the (r, z) plane showing simulation R\I3 - CHI20 at 
t = 5000 yr. The black/blue contours enclose the computational cells con- 
sisting of the original ejecta/ shrapnel material by more than 90%. The color 
bar indicates the logaritmic density scale and the linear temperature scale. 
Right panel: same as left panel, fort= 1 1000 yr. 

Rayleigh-Taylor and Richtmyer-Meshkov instabilities are visible 
as finger-like structures both in the density and temperature maps. 
At this stage, the knot is partially eroded by the hydrodynamic in- 
stabilities and evolves toward a core-plume structure. The core of 
the knot, however, is still significantly overdense with respect to 
the surrounding shocked ejecta. The right panel of Fig. |2] shows 
the shrapnel at ^ = 11000 yr, with its characteristic supersonic bow 
shock protruding beyond the SNR main shock. 

WC02 found that ejecta knots with density contrast;^ < 100 
are rapidly fragmented and decelerated in the inter shock region and 
do not even reach the main shock front (these eff'ects being more 
dramatic for small clumps). Nevertheless, we notice that the value 
of X in WC02 refers to the onset of the interaction between the 
knot and the reverse shock and that x is not constant during the 
evolution of the system. In the "free" expansion phase, the den- 
sity of the shrapnel does not drop down uniformly (as that of the 
other ejecta does) and the shrapnel undergoes both diff'usion and 
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-1 

r (10^® cm) 



Figure 3. Density 2-D cross-sections through the (r, z) plane for model 
Rl/3 - CHI20 at t = 2500 yr. The color scale increases linearly between 
10~^^ g/cm^ and 2 x 10~^^ g/cm^. The contours enclose the computational 
cells consisting of the original shrapnel material by more than 99% (bold 
line) and 50% (narrow line). 
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Figure 4. Temporal evolution of the shrapnel /ejecta density contrast for 
model Rl/3- CHI20. The blue horizonthal line marks the beginning of the 
interaction between the shrapnel and the SNR reverse shock 



expansion. Figure [3] presents a close-up view of the shrapnel den- 
sity structure at f = 2500 yr, showing that, while the outer parts 
of the knot diffuse and mix with the expanding ejecta, its central 
core remains much denser. The density of the core of the clump 
drops down much more slowly than that of the spherically expand- 
ing ejecta. Therefore, the inhomogeneous rarefaction of the knot 
makes the density contrast between the core of the shrapnel and the 
expanding ejecta higher, and rapidly increases until the shrapnel 
reaches the reverse shock. We computed x during the expansion 
phases, by calculating the shrapnel density, p7, as the average of 
the density in all the computational cells where the shrapnel con- 
tent is > 909^1. We then divided ps by the ejecta density (along the 
r axis) at the same distance from the origin as the shrapnel cen- 

^ These cells are closer to the center of the knot, and are less affected by 
the diffusion. 



ter. Figure U] shows the evolution of ;^ as a function of time for the 
Rl/3 - CHI20 simulation. The figure shows that the ejecta knot 
reaches x > 100 as it approaches the reverse shock, hence our re- 
sults are in agreement with those of WC02. Our model shows that 
a knot that was only 20 times denser than the surrounding ejecta 
(at the beginnig of the simulations) can reach the SNR main shock 
without being fragmented in the intershock region and can produce 
protrusions that are similar to those actually observed in the Vela 
SNR. 



3.2 Effects of thermal conduction and radiative cooling 

Figure |5] shows the 2-D cross-sections through the (r, z) plane of 
temperature and density at ^ = 5000 yr and t = 11000 yr for the 
R\/3- CHI20 - TR simulation (same parameters as /?1 /3 - C///20, 
but including radiative cooling and thermal conduction). The diffu- 
sive thermal conduction completely suppresses the formation and 
the development of hydrodynamic instabilities and smoothes the 
temperature and density profiles. This result is in agreement with 
expectations, as shown below. The characteristic amplitude growth 
rate, da/dt of a single-mode perturb ation of Richtmy er-Meshkov 



instabilities can be calculated as (see 

^ - kAvaA 
dt 



) ation 01 Kicntmy e 
lRichtmvei]ll96Ql) 



(1) 



where k is the perturbation wavenumber, Av is the velocity jump at 
the instability and A - {pi - p2)/(p\ + Pi) is the Atwood number. 



The characteristic time- scale, 
the perturbation is therefore 

/ 



-= [s] 



IttAvA 

where / is th e structure size. 
ISpitzeilll962h . 



a /(da/dt), for the growth of 



(2) 

As for the thermal conduction (see 



dE\ 2 T 

— ^V-[k{T)VT]- -K{T)- 

Icond 7 P 



(3) 



where k(T) = 5.6 x 10~^r^/^ erg s~^ cm~^ is the Spitzer's co- 
efficient and / is the characteristic length of temperature variation, 
therefore, the thermal conduction time- scale is: 



Ink f o ^ 1 
^-''=5(731)4?) -2.6x10 ^ 



[s] 



(4) 



For a characteristic structure with size / ~ 4 X 10^^ cm, parti- 
cle density n - 0.% cm"^, Atwood number A ~ 0.44, Av ~ 2 x 10^ 
cm/s, T ~ 1.3 X 10^ K (similar to that shown in Fig. [J]), Tinst ~ 
^■^Tcond ~ 2000 yr. Therefore, the thermal conduction diffusive 
processes develop faster than the hydrodynamic instabilities and 
density and temperature inhomogenieties are smoothed out before 
they can grow. 

The evolution of the position of the shrapnel head and the pro- 
trusion it produces to the remnant shock front are similar to those 
obtained without including thermal conduction and radiative cool- 
ing. Nevertheless, the shrapnel evolution is remarkably different 
from that obtained in the pure HD simulations. In particular, as 
shown by the blue contours in Fig. \5\ the ejecta knot is elongated 
along its direction of motion and rapidly assumes a cometary shape, 
characterized by a prominent tail which is rich in shrapnel material. 
After t = 11000 yr, shrapnel material is present at ~ 10 pc away 
from the shrapnel head. Moreover, the shrapnel material is effi- 
ciently heated by thermal conduction with the surrounding shocked 
ejecta. Let MXshra be the mass of the plasma in the computational 
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Figure 5. Same as Fig.|2]for model Rl/3 - CHI20 - TR. 



cells consisting of the original shrapnel material by more than 90% 
and having a temperature higher than 10^ K (and therefore emit- 
ting thermal X-rays). Figure |6] shows the evolution of MXshra as 
a function of time for the simulations i?l/3 - CHI20 - TR and 
R\/3 - CHI20. When thermal conduction is at work, ~ 90% of 
the original shrapnel mass is heated up to X-ray emitting temper- 
ature at t = 11000 yr (i. e., at the age of the Vela SNR), while 
if thermal conduction in inhibited, only 50% of the original mass 
has temperature higher than 10^ K. Figure[6]also shows the amount 
of X-ray mass beyond the shock front. In the pure HD simulation 
the hot shrapnel material is all beyond the SNR shock front. In the 
Rl/3 - CHI20 - TR simulation, only part of the X-ray emitting 
shrapnel is beyond the shock front and there is a significant frac- 
tion of the ejecta knot material (the shrapnel tail) that is inside the 
SNR shell and that is expected to emit thermal X-rays (see Sect.|4l). 

We point out that in SNRs the efficiency of thermal conduc- 
tion can be significantly reduced by the presence of the magnetic 
field (which is not taken into account in our model). If we as- 
sume an organized ambient magnetic field, the thermal conduc- 
tion is anisotropic, because the conductive coefficient in the direc- 
tion perpendicular to the field lines is several orders of magnitude 
lower than that parallel to the field lines, which coincides with the 
Spitzer's coefficient k(T). The effects of the magnetic-field-oriented 
thermal conduction in the interaction betw een shocks and dens e 
clump have been investigated in detail in lOrlando et al.l (l2008 l). 
Because of the high beta of the plasma, the magnetic field lines 
are expected to envelope the hydrodynamic fingers thus hampering 
the thermal conduction with the surrounding material. At the same 
time, the magnetic field is expected to be trapped at the top of the 
ejecta clumps, and this yields to an increase of the magnetic pres- 
sure and field tension w hich limits the gr owth of hydrodynamic in- 
stabilities (see 012 and lSano et al.ll2012h . The basic physics of the 
interaction between the ejecta knots and the SNR shocks is simi- 
lar to that for the interaction of planar shocks with an interstellar 
cloud (e. g. WC02) and it has been shown that, in this case, sim- 
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Figure 6. Temporal evolution of MXshra (mass of the plasma in the com- 
putational cells consisting of the original shrapnel material by more than 
90% and having a temperature higher than 10^ K) for the simulations 
Rl/3 - CHm - TR (blue curves) and Rl/3 - CHm (black curves). The 
dashed lines indicate the part of MXshra that protrudes beyond the SNR 
shock front. 



ulations including thermal conduction in an unmagnetized plasma 
and pure HD simulations are limiting cases that encompass the re- 
sults obtained with diff'erent configurations of the magnetic field 
dOrlando et aLl l2008Y We can then conclude that our simulations 
provide the two extreme cases that bracket all the possible interme- 
diate scenarios. 



3.3 Effects of the initial conditions 

We study the eff'ects of the initial conditions on the shrapnel evolu- 
tion with different simulations, as shown in Table [T] In particular, 
we explored two diff'erent initial positions of the shrapnel in the 
ejecta profile (R 
contrasts (x = 10, 20, 50) 

In agreement with WC02, we found that shrapnel formed in 
the inner ejecta layers (i. e., those that reach the reverse shock later) 
produce smaller protrusions. In particular, an ejecta knot originat- 
ing at Rs = 1/6 R^^j, does not even reach the SNR shock in the 
time spanned by our simulations. Therefore, our model indicates 
that shrapnel A-F (all protruding well beyond the Vela main shock) 
originated in more external layers. Left panel of Fig. |7] shows the 
2-D cross-sections through the (r, z) plane of temperature and den- 
sity at the age of the Vela SNR for the Rl/6- CHI50 simulation. In 
this case, the knot is well within the intershock region, even though 
its initial density contrast (x = 50), was higher than that of the 
Rl/3 - CHI20 run. However, our models of knots originating at 



1/6 R^^j, 1/3 R^^j) and three different density 



Rs = 1/3 R^^j clearly show that denser shrapnel produce deeper 
protrusions and are more stable against the fragmentation and the 
deceleration induced by the hvdrodynamic instabilities in the inter- 
shock region (as in WC02jj. Fig. [T] right panel, shows the 2-D 
cross-sections through the (r, z) plane of temperature and density 
at the age of the Vela SNR for the Rl/3- CHI50 simulation. As ex- 
pected, the shrapnel head is much further away from the SNR shell 



^ In our runs, higher values of correspond to smaller values of the shrap- 
nel cross-section (given that the shrapnel mass is fixed to 1.19 X 10^3 

g in 

all our simulations), and this concurs in making the knot a more penetrating 
bullet. 
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-24.75 -22.2S, 
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Figure 7. Left panel: density (left) and temperature {right) 2-D cross- 
sections through the (r, z) plane showing simulation Rl/6 - CHI50 at 
t = 11000 yr. The black/blue contours enclose the computational cells con- 
sisting of the original ejecta/ shrapnel material by more than 90%. The color 
bar indicates the logaritmic density scale and the linear temperature scale. 
Right panel: same as left panel, for simulation Rl/3 - CHI50. 
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Figure 8. Same as Fig.|4]for model Rl/3- CHI50. 



than in the R\/3- CHI20 case (shown in the right panel of Fig.|2]). 
We computed the time evoulution of ;^ for the Rl/3 - CHI50 run 
by following the procedure described in Sect. 13.11 We found that, 
in this case, the ejecta knot reaches the reverse shock with a very 
high density contrast (x > 1000, see Fig. [8]), thus producing a very 
prominent protrusion in the SNR. 

Figure |9] shows the position of the shrapnel head (in units 
of the shell radius) as a function of time for all our simula- 
tions. Figure [9] also sho\y s the projected positions of Shrapnel A-D 
(lAschenbach et al.ll 19951) and of the ejecta knots FilE and RegNE 





aO 


- 










^ ' RegNeJi 




$FilE - 










D 2 4 e 


8 10 12 



Figure 9. Position of the shrapnel head (in units of the shell radius Rsheii) as 
a function of time for all the simulations listed in Table[T] Black curves indi- 
cate pure hydrodynamic models, while blue curves indicate simulations in- 
cluding thermal conduction and radiative cooling. Dashed/solid/thick lines 
indicate;^ = 10/30/50, respectively. The red diamonds show the projected 
positions of Shra pnel A-D ( Aschenbach et al. 1995) and of the ejecta knots 



positions 01 Shra pnel A-D ( Aschenba 
FilE and RegNE iMiceH et al.ll2008h . 



("Miceli et al'lOOS") with respect to the position of the shock front in 
the Vela SNR. These values were calculated by approximating the 
Vela SNR as a circular shell with angular radius 211' and center 
with coordinates aj2ooo = 8^^36'"19.8^ 6j2ooo = -45°24'45". Mod- 
els including the effects of radiative cooling and thermal conduction 
(blue curves in Fig. ^ do not provide significant diff'erences with 
respect to pure HD models (black curves) in terms of the shrapnel 
position. 

By considering all the simulations with Rs = 1/3 R^^j, we find 
that the position of the head of the knot, Rh, after ~ 1 1000 yr ranges 
between ~ 1.1 Rsheii (ior x =10) and ~ 1.4 Rsheii (ior x - 50). 
These values are similar to those observed for Shrapnel B, C, and 
D. 

Shrapnel A, E, and F (Shrapnel E, F are not shown in Fig. |9) 
instead, have = 1.51 Rsheih Rh = 1.88 7?^/^^//, and /?/^ = 1.91 Rsheih 
respectively. These values are much larger than those obtained in 
our simulations. Our results clearly suggest that these large dis- 
tances from the shell can be produced by ejecta knots originating 
in outer ejecta layers. An alternative possibility is that the original 
density contrast for these shrapnel was > 50. Nevertheless, Fig. |9] 
shows that the final position of an ejecta clump is more sensitive to 
its initial position and depends only weakly on x- In fact, by vary- 
ing the initial position of the knot by a factor of two, we found that 
its final position varies approximately by the same factor, while a 
variation of the initial density contrast by a factor of five, only de- 
termines a 30% variation in the final position. Therefore it is more 
likely that Shrapnel A, E, and F were produced at 7?^ > 1/3 R^^j. 

The two simulations with Rs = 1/6 R^^j show that the ejecta 
knots originating in the inner ejecta layers do not reach the forward 
shock, even for the highest density contrast (x = 50). The position 
of the head of the knot predicted b y our simulation s is in agree- 
ment with that observed for FilE (Miceli et al ] |2008 h. X-ray emit- 
ting ejecta knots have been obser ved, in projection, inside the Vela 
SNR shell (e. g., RegNE and FilE jMiceli et alJ2008l) . However, we 
point out that the actual position of these "internal" shrapnel can be 
well outside the SNR shell, and therefore the values reported in Fig. 
[9] might be considered as lower limits. 
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4 DISCUSSIONS AND CONCLUSIONS 

The structure of the ejecta in a SNR contains the imprint of the 
metal-rich layers inside the progenitor star, and may help to under- 
stand the processes occurring in the latest stage of stellar evolution. 
We performed a set of hydrodynamic simulations to study the evo- 
lution of the ejecta knots in the Vela SNR. 

We found that moderately overdense clumps (initial density 
contrast;^ ~ 10) can produce protrusions in the SNR shell simi- 
lar to those observed for the Vela shrapnel. WC02 found that only 
clumps that reach the reverse shock with density contrast ;^ ^ 100 
can reach the main shock front and produce significant protrusions. 
This criterium is fulfilled in all our simulations. In fact, the (ini- 
tially) moderately overdense clump experiences diff'usion in the 
"free" expansion phase, and, while its outer parts mix with the 
surrounding ejecta, its central core remains much denser. This in- 
homogeneous rarefaction makes the density contrast between the 
core of the shrapnel and the expanding ejecta higher as the remnant 
evolves, and;^ reaches values ~ 100 - 1000 at the interaction with 
the reverse shock. 

In particular, a knot with initial;^ = 20 and Rs = 1/3 R^^. (sim- 
ulations R\l?> - CHI20 and R\/3 - CHI20 - TR) can explain the 
observed features associated with the Vela Shrapnel D. Figure [TOl 
shows the the ROSAT All-Sky Survey image of the Vela Shrap- 
nel D in the 0. 1 - 2.4 keV energy band, compared with a synthesis 
of the X-ray emission in the 0.1 - 2.4 keV band derived from the 
Rl/3- CHI20 - TR simulati on at 11000 yr. The synthesized X-ray 
map has been obtained as in lMiceli et al.l (l2006h : we produced the 
3D map of the emission measure and temperature in the cartesian 
space (x', y, z'), where the y' axis corresponds to the direction of 
the line of sight and is perpendicular to the (r, z) plane. We de- 
rived the distribution EM{x\ z') vs. T{x\ z') by considering all 
the contrbutions along the line of sight, for each {x' , z'). We then 
synthesiz ed the map of the X-ray emission using the MEKAL spec - 
tral code (iMewe et al.ll985LlMewe et al.ll986l . lLiedahl et al.ll995h . 
assuming a distance of 250 pc, and an interstellar column density 
Nh = 2 X 10^*^ cm~^. Finally, we degraded the spatial resolution of 
the synthesized X-ray map to match the resolution of the ROSAT 
image and randomized the map by assuming Poisson statistics for 
the counts in each image bin. 

Figure [TO] shows that the bright X-ray spot at the apex of the 
shrapnel bow- shock and the enhanced X-ray luminosity at the base 
of the protrusion can be explained by our model as local enhance- 
ments of the plasma emission measure. The overall shape of the 
observed features is also reproducted by our model. Moreover, the 
mass of the X-ray emitting shrapnel predicted by our model is in 
agreement with that measured for Shrapnel D. In fact, it has been 
estimated that the X-ray emitting mass (consider ing only the bulge 
above the Vela SNR bow shock) is ~ O.IM© OKatsuda & Tsunemj 
l2005h i. e., the same order of magnitude as that obtained in our 
simulation (see dashed curves in Fig. [6]). Vela Shrapnel D is there- 
fore consistent with being originated by an ejecta knot 20 times 
denser than the surrounding ejecta and initially located at 1/3 R^^y 
As shown in Fig. |9] we found a degeneracy in the space of the 
initial parameters, and similar results can be obtained also by en- 
hancing/decreasing the density contrast and decreasing/increasing 
the initial distance of the knot to the explosion center. However, as 
explained in Section [331 the final position of the shrapnel is much 
more sentitive to its initial position in the ejecta profile. In conclu- 
sion, according to our model, the original position of Shrapnel D 
should not diff'er much from 1/3 R^^.. 

As shown in Sect. 13.31 our model suggests that Shrapnel B, C, 




Figure 10. Upper panel: Rosat All-Sky Survey image of Vela Shrapnel D in 
the 0.1 - 2.4 keV energy band. The bin-size is 1.5' and the image has been 
smoothed through a convolution with a Gaussian with cr = 3 pixels. North 
is up and East is to the left. Lower panel: Synthesized X-ray emission in the 
0.1 - 2.4 keV band for the simulation /?l/3 - CHm -TRatt= 11000 yr 
(see Sect.|4j. 
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and D were all originated at ~ 1/3 R^^.. This conclusion is in agree- 
ment with the results of X-ray data analysis t hat show that Shrap- 
nel B and D have s imilar abundance patterns ( Katsud a & Tsuneniiil 
l200llYamaguchi & Katsuda 2009). The fragment RegNE, which 
appears inside the V ela shell, has similar abundances as Shrapnel D 
(iMiceli et"aDl2008h and its position is compatible with an origin in 
the same layer where Shrapnel B, C, and D were generatecJl. 

In Sect. 13.31 we also pointed out that it is highly unlikely 
that Shrapnel A originated in the same ejecta layer as Shrap- 
nel D. Indeed, the abundance patterns observed in Shrapnel A 
are remarkably different fr om that observed in Shrapnel B and D 
(iKatsuda & Tsunemi|[200^ . thus suggesting a different location of 
this knot in the ejecta profile. Nevertheless, the Si:0 ratio is much 
higher in Shrapnel A than in Shrapnel D and this indicates that 
Shrapnel A comes from a region of the progenitor star below that 
of the Shrapnel D (e. g. iTsunemi & Katsudd|200 6^). at odds with 
our predictions. As explained above, an unrealistically high initial 
X is required for an inner shrapnel to overcome outer knots, if we 
assume that the initial velocity profile of the ejecta increases lin- 
early with their distance from the center. A possible solution for 
this puzzling result is that the Si burning layer (or part of it) has 
been ejected with a higher initial velocity, e. g., as a collimated jet. 
It is noteworthy to remark that in other core-collapse SNRs the Si- 
rich ejecta may show a very peculiar jet-counterjet structure. The 
well known case of Cas A has been studied in detail thanks to a very 
long Chandra observation (Hwang et al. 2004) showing a jet (with 
a weaker counter] et structure) composed mainly of Si-rich plasma. 
iLaming et al.l (l2006h have performed X-ray spectral analysis of sev- 
eral knots in the jet and concluded that the origin of this interesting 
morphology is due to an explosive jet and it is not arising because 
of an interaction with a cavity or other ISM/CSM peculiar structure. 
Therefore a jet origin for the Si-rich knots is sound. 

Finally, we investigated the effects of thermal conduction, 
finding that it determines an efficient "evaporation" of the ejecta 
knot and accelerate its mixing with the surrounding medium. More- 
over, it affects the shrapnel morphology, producing the formation of 
a long, metal-rich tail. Enhanced metallicities have been observed 
in the tails of Shrapnel A, B, and D and the abundance analysis per- 
formed on the X-ray spectra clearly suggests an efficient mixing of 
the ejecta knots with the surrounding medium (Katsuda & Tsunenii! 
[2OO5IJ 2OO6. and Yamag uchi & Katsuddl2009 ). These results are in 
qualitative agreement with our findings. A quantitative compari- 
son between models and observations requires a forward model- 
ing approach, consisting in the synthesis of the X-ray spectra from 
the simulations and a detailed comparison between the synthesized 
observables and the corresponding obse rvations (e. g., thr ough a 
spatially resolved spectral analysis, as in iMiceli et al.ll200a . e. g.). 
Moreover it will also be important for future models to include 
some seed magnetic fields in the simulations to study how tangled 
the field becomes and what this implies for thermal conduction. 
These further studies are beyond the scope of this paper. 

In conclusion, our hydrodynamic modelling of ejecta knots 
in the Vela SNR allowed us to find that: i) the observed shrapnel 
can be the results of moderate density inhomogeneities in the early 
ejecta profile; ii) the evolution of a shrapnel in the SNR is very sen- 
sitive to its initial position and depends much less (but does depend) 
on the initial density contrast; iii) thermal conduction plays an im- 



^ In this case, its projected distance from the center of the SNR must be 
much smaller than its actual distance. 



portant role and explains the efficient mixing of the ejecta knots 
observed in X-rays. 
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APPENDIX A: TEST ON SPATIAL RESOLUTION 

As explained in Sect. 12.11 the spatial resolution of our simulations 
is reduced by a factor of 2 at ? > 2500 yr. To check whether this res- 
olution is sufficient to capture the basic evolution of the system, we 
repeated simulation Rl/3 - CHI20 by mantaining the finest level 
of resolution at its initial value (1.95 X 10^^ cm, hereafter simula- 
tion Rl/3 - CHI20 - hiRES) for the whole run. We verified that 
simulations Rl/3 - CHI20 and Rl/3 - CHI20 - hiRES yield very 
similar results and all the quantities discussed in the paper (posi- 
tion of the shrapnel as a function of time. X-ray emitting mass of 
the shrapnel, etc.) do not change significantly. Pigure lAll shows 
the 2-D cross-sections of the density through the (r, z) plane ob- 
tained for simulations Rl/3 - CHI20 and Rl/3 - CHI20 - hiRES 
at f = 5000 yr. This evolutionary stage is the most critical, since 
the system is still relatively small and the ejecta knots interact with 
the hydrodynamic instabilities in the intershock region. Pigure lAll 
clearly prove that simulation R1/3-CHI20 already provides a very 
accurate description of the system and run Rl/3- CHI20 - hiRES 
does not introduce any major diff'erences. 



